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o : 

f - *^ Abstract. We study the short time heat content asymptotics for spectral 

, boundary conditions. The heat content coefficients are shown to be non-local 

and some preliminary results concerning the structure of the first few terms 
are given. 

<D ■ 

m : 
\o ■ 

1. Introduction 

t-H . 

^ • The heat trace asymptotics for 'exotic' boundary conditions have recently at- 

| tracted considerable interest; for a survey over this field see [Q. This article is 

devoted to the analogous questions for the heat content asymptotics. Whereas 
some formulas are available for the heat trace asymptotics with spectral boundary 
conditions ||, [l0| , nothing is known about the heat content asymptotics in this 
CsJ | setting. We shall postulate in equation (J2J) the existence of an appropriate asymp- 

totic series and then use special case calculations of spinors on the unit ball to study 
C"| \ the heat content asymptotics defined by non-local spectral boundary conditions - 

Cln' these are preliminary results in an ongoing investigation. 

Spectral boundary conditions were first introduced in the study of the index the- 
orem for manifolds with boundary by Atiyah, Patodi, and Singer [|l] who assumed 

S that the structures involved were product near the boundary. Their work was later 
extended by Grubb and Seeley |1| |l6], |l7], |l8|, [l9| to the general setting, 
^s! I We briefly establish the notational conventions we shall employ and refer to [|[ |lo| 

for further details. Let M be a compact m-dimensional Riemannian manifold with 
smooth boundary clM. We suppose given unitary vector bundles Ei over M. and 
■ a first order elliptic complex: 

P : C°°{Ex) -> C°°(E 2 ). 

As such an elliptic complex need not admit local boundary conditions, it is nat- 
ural to impose non-local spectral boundary conditions, which may be described as 
follows. Let 7 be the leading symbol of P. Let V be an auxiliary unitary con- 
nection with V7 = 0; in many applications there is a natural choice available, but 
it is convenient to work quite generally for the moment. Let V m denote covariant 
differentiation with respect to the inward geodesic unit normal. Near the boundary, 
we decompose 

P=7m{V m + S}, 
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where B is the associated tangential first order operator on C°°(-Ea lax). We also 
suppose given an auxiliary self-adjoint endomorphism of E\\qm which we use to 
define a self-adjoint tangential operator A on C co {Ei\qm) by setting: 

A= ±{P + P*} + 0; 

here the adjoint of B is taken with respect to the structures on the boundary. The 
endomorphism arises naturally in the study of the signature and spin complexes 
where the metric is not product near the boundary and in this context is expressible 
in terms of the second fundamental form ; compensates for non-canonical choice 
of V m and plays an important role in the analysis in Section |3|. 

Let II be orthogonal projection on the span of the eigenspaces corresponding 
to the non-negative eigenvalues of A. Let Pn be the realization of P with the 
boundary condition II. The index theorem for manifolds with boundary expresses 
index(Pn) in terms of characteristic forms integrated over M., a compensating in- 
tegral involving the second fundamental form over dM. , and a non-local term (the 
eta invariant) for the classic elliptic complexes [0, 

To simplify the discussion, we shall assume that E\ = E2 = E, that P is formally 
self-adjoint, and that P is of Dirac type; thus D := P 2 is formally self-adjoint and 
of Laplace type. The leading symbol 7 of P is skew-adjoint and 

7 (£) 2 = -|£| 2 Id for any £ e T*M. 

Db is the realization of D defined by the operator 

BV' := n^| aA1 © ILP</W 

To avoid some technical fuss with the zero mode spectrum, we shall suppose that 
ker{j4} = {0}. We shall also assume that 

(1) 1mA = -A-f m so 7„ l II = {Id - II}7 m . 

The operators Pn and Dg = P^ are then self-adjoint. 

We now describe the fundamental solution of the heat equation. Let /1 rep- 
resent the initial "temperature" distribution of the manifold M.. The subsequent 
temperature distribution h(x,t) for t > is then given as the unique solution of 
the equations: 

dth(x,t) = Dh(x,t), Bh(x,t) = Q, and lim h(x,t) = fi(x). 

Let (•, •) denote the Hermitian inner product on E. Let dx and dy be the Riemann- 
ian measures on M. and dA4, respectively. Let fa represent the "specific heat" of 
the manifold. The total heat energy content of the manifold is given by: 

0{hJ 2 ,D,B) (t) = J M dx (h(x, t) , h (z)> ■ 

It is convenient at this point to review the situation for 'standard' boundary 
conditions. Let D be a formally self-adjoint second order operator of Laplace type 
whose realization is defined by Dirichlet or Robin boundary conditions B. It is 
known that there is a complete short time asymptotic series of the form: 

00 

P(h,f 2 ,D, B) ~ Yl Wuh, D, B)t n ' 2 . 

71=0 

The coefficients f3 n are called the heat content asymptotics. 

Let Roman indices {i, j, fc, 1} (resp. {a, b}) range from 1 to to (resp. to — 1) and 
index an orthonormal frame for the tangent bundle TM. (resp. TdM). Let L a i> be 
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the components of the second fundamental form and let Rijki be the components of 
the Riemann curvature tensor (with the sign convention Ri 22 i = +1 on S 2 C R 3 ). 
Let ';' be multiple covariant differentiation with respect to the Levi-Civita connec- 
tion on M. and the natural connection defined by D, and let £ be the endomorphism 
of E defined by D 01. We adopt the Einstein convention and sum over repeated 
indices. We refer to ||] for the proof of the following Lemma: 

Lemma 1. Let D be an operator of Laplace type on a compact Riemannian mani- 
fold M. with smooth boundary dM. . 

(1) If Btp = ip\dM defines Dirichlet boundary conditions, then 

(a) l3 Q (fi,f 2 ,D,B) = J M dx (f u f 2 ). 

(b) [3 1 (f 1 ,f 2 ,D,B) = -2Tr- 1 / 2 f dM dy (hJ 2 ). 

(c) (3 2 (fi,f2,D,B) = -f M dx (Dfi,h) 

+ JdM d V {(l L aafl,f 2 ) - (fl, f 2 - m )} 

(d) Ps(fi,f2, D, B) = -27T-V3 J M dy {§</i ;mm , f 2 ) + / 2;mm ) 

— {fl\a,h;a) + (£f,f 2 ) — § L aa (f\- m , f 2 ) — ^L aa (f 1 , f 2 . m ) 
-^((j^LaaLbb ~ gLabLab + g Ramam)fl , f 2 ) } ■ 

(2) If Btp = (V m + S)ip\gM defines Robin boundary conditions, then: 

(a) p (h,f2,D,B)=J M dx(f 1 ,f 2 ). 

(b) p 1 (f 1 ,f 2 ,D,B) = 0. 

(c) l3 2 {h,f 2 ,D,B) = -$ M dx (Df l7 f 2 ) + f dM dy {BfrJ 2 ). 

(d) [3 3 (f 1 ,f 2 ,D,B) = l-ir~ 1 / 2 f dM dy (Bfx,Bf 2 ). 

Lemma |l| shows that the heat content coefficients [3 n for n < 3 are locally 
computable for Dirichlet and Robin boundary conditions. In fact, all the coefficients 
[3 n are locally computable for these boundary conditions. These invariants have 
been studied extensively [||, f| ||, ||^, ; we refer to (l| for a recent survey article. 

We now return to the setting of spectral boundary conditions. We shall assume 
the existence of a similar asymptotic series for the realization of D = P 2 defined 
by spectral boundary conditions as t — > 0: 

oo 

(2) P(fiJ 2 , D,B)~J2 Wu h,D, B)t n / 2 . 

71 = 

A natural question then arises, why do we believe there are no log terms - after 
all, in the expansion of the heat trace Tr L 2(e~ tL>B ), log terms appear! Our answer 
is two-fold. First of all, the heat content asymptotics exhibit various structural 
simplifications compared to the heat trace asymptotics. One example is, as we 
shall see presently in Section Ej, that in the special case of the Dirac operator on 
the unit ball the relevant universal constants do not depend on the dimension, as 
they did for the heat trace aymptotics. Another example is, that the heat content 
asymptotics do not show any signs of the loss of strong ellipticity for the case of 
oblique boundary conditions [JlTJ . Furthermore, the heat trace asymptotics give 
rise to log terms above the dimension of the manifold in the series; for n < m, the 
asymptotic coefficients do not have log terms. So even if the situation is as for the 
heat trace, our results will still hold true for n < m. But this is certainly a question 
that merits further investigation. 

The first coefficient (3q is easily described. Since limt_>o h(x, t) = fi(x), we have 
limt^o /3(/i> Hi D,B)(t) — f M dx (fi,f 2 ) and thus, as for Dirichlet and Neumann 
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boundary conditions, 

(3) Wuh,D,B) = f M da(fi,h). 

Here is a brief outline to this paper. In Section^, we discuss functorial properties 
of these invariants, show they are non-local, and outline what we believe the formula 
for /3\ and (3 2 to be. These results are based on the special case computations in 
Section || giving a complete calculation of the heat content function for the Dirac 
operator on the unit ball with the standard metric. 

2. Functorial Properties 

The invariants (3 n for Dirichlet and Robin boundary conditions have a number 
of functorial properties Q which extend immediately to this setting: 

Lemma 2. Let = P^ be a self-adjoint operator of Laplace type defined by 
spectral boundary conditions as described above. Then: 

(1) We have (3 n {fu h,cT 2 D, B) = c™- n f3 n {h, / 2 , D, B) for any c e M+. 

(2) We have /3„(/ 1; f 2 ,D,B) = /3„(/ 2 , fi,D, B). 

(3) 7/6/1 = 0, then[i n {f l J 2 ,D,B) = -ll3 n -i{Df l ,h,D,B). 

Proof. The operator Dq has a discrete spectral resolution Sd,b '■= {ip u ,\v} [ [l5| 
with associated Fourier coefficients: a v (tp) := f M dx (ip,ip v ). We may then express: 

(4) h{x,t)=Y, v e- tK °u{h)i>»{x) so 

0(fx, h, D,B)(t) = E„ e-^a v {h)au{h). 

The estimates of show these series converge uniformly. 

We take Q(c) := c -1 to ensure the associated boundary condition is unchanged. 
Note that the Riemannian measure defined by the operator c~ 2 D is c m dx. Since 
S c -*d,b = {c-™/ 2 7/v,c- 2 A„} and <(/) = c m / 2 a„(/), we have: 

fr, c~ 2 D, B)(t) = c m (3(f u f 2 , D, B){c~ 2 t). 

Assertion (1) follows by equating powers of t in this equation. Since the roles of 
fi and / 2 are symmetric in display (|]), assertion (2) follows. If Bf\ = 0, then the 
boundary terms vanish and we can integrate by parts to compute: 

(?v(Dfi) = J M dx {Dfx^v) = f M dx {fi,D%l) y ) = \„a v (fi) so 
-d t /3(hJ 2 ,D,B)(t) =Y,»^- tx »o- v {h)o- v {h) = f3{Dh,f 2 ,D,B){t). 

We can now establish (3) by equating terms in the asymptotic expansions for 
dtP(fij2,D,B)andp(Df 1 ,f 2 ,D,B). ' □ 

The following is an important observation. 

Lemma 3. The heat content coefficients for spectral boundary conditions are not 
locally computable 

Proof. If Bfi = then Px{f\, / 2 , D, B) = by Lemma ||. If f3i is locally computable, 
then dimensional analysis (i.e. the scaling property given by assertion (1) of Lemma 
implies that there must exist a universal constant Co(ro) so that: 

/?i(/i> h, D, B) = 2-k- 1 ' 2 J dM dy c (m) (/i,/ 2 ). 

Since generically there are eigensections with Bfi = but fi\dM ^= 0, we must have 
c = 0; so far, the argument is exactly the same as for Robin boundary conditions 
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given in Q to prove Lemma [l] (2b). However, the calculation on the ball that 
we shall present in Section |^ shows the power t 1 / 2 is present in the asymptotic 
expansion with spectral boundary conditions. This contradiction establishes the 
Lemma. □ 

The coefficient (3q is given by equation (Q). Using the principle of not feeling 
the boundary, we see that the interior integrals defining [3 n for spectral boundary 
conditions are the same as those defining (3 n for either Dirichlct or Robin boundary 
conditions. Writing the interior term for (3 2 in the form (-D/1,/2) destroys the 
symmetry of Lemma || (2) so instead we use (P/i, P/2) and add suitable boundary 
correction terms. The boundary operator II is a 0-th order operator; it is unaffected 
by rescaling. To ensure that properties (1) and (3) of Lemma |^ are satisfied, i.e. 

f3 n (hJ 2 ,c- 2 D,B) = c m - n (3 n (fuf2,D,B) and 

Wi,f2,D,B)=0 n (h,h,D,B) 

we are lead to consider the following ansatz for [3\ and /3 2 - the only non-local terms 
are introduced by the boundary condition. 

Ansatz 4. There exist universal constants Ci = Ci(m) so that: 

(1) f3 1 (f 1) f 2 ,D,B) =2tt- 1 / 2 J gM dy c (m) (n/ lf n/ 2 ). 

(2) f3 2 (h,f 2 , D,B) = - J M dx (PA, Pf 2 ) + J dM dy { Cl (m)(n 7m PA, Uf 2 ) 

+ Cl (m)(n/ 1 ,n 7m p/ 2 ) + C2 (m)L M (n/ ll n/ 2 ) + C3 (m)(en/ 1 ,n/ 2 )}. 

The remainder of this note is devoted to the evaluation of these coefficients: 

Lemma 5. We have co(m) — —1, Ci(m) = 1, c 2 {m) — ^, and €3(171) = 0. 

It is interesting that these coefficients seem to be dimension free; the usual 
trick of dimension shifting employed in || does not work with spectral boundary 
conditions. By contrast, the coefficients in the heat trace asymptotics for spectral 
boundary conditions are highly dimension dependent || [HJ . 

Proof. We must ensure the properties of Lemma 2 are satisfied. In particular 
/M/i, /2) D, B) must vanish if II/ 1 =0. The boundary term for 0x must be homo- 
geneous of degree 0. Since II is a th order operator, (n/i,n/ 2 ) has the correct 
homogeneity, is symmetric in the roles of {/i,/2}, and vanishes when Bf\ = 0. 
This motivates the formula given in (1). 

Consider j3 2 (fi, f 2 ,D,B). To preserve the interior symmetry, we use the interior 
integrand — (P/i, P/2) rather than —(D/1,/2). The corresponding integrals are 
related by the formula: 

- J M dx{Dh,f 2 ) = - J M dx{Pf u Pf 2 ) + J dM d y ( lrn Ph, f 2 ). 
If Bfi = 0, we know from Lemma || (3) that 

ft(/iJa.AB)=-A (Dfi ,f 2 ,D,B). 
Since II/i = ILP/i = on dM, we use equation (]l|) to see that 
ImPfx - 7m (1 - H)P/i = n 7m P/x on dM. 
Consequently, we find: 

- ftP/i,/ 2 , D,B) = - j M dx (P/i, P/ 2 ) + J dM dy (n 7m P/!, f 2 ) 

(5) = -S M dx{Ph,pf 2 ) + S dM dy{n lm Ph,nf 2 ). 



6 



PETER GILKEY 1 , KLAUS KIRSTEN 2 , AND JEONG-HYEONG PARK 3 



This shows ci(m) = 1. Lemma || (2) shows we need to include ci(m)(Hfi, Hj m Pf2) 
into Ansatz ^| (2). Lemma || (3) shows that apart from the invariants multiplied by 
ci(m) additional invariants must disappear if LT/i = HPfi = 0. Applying Lemma 
H (1) this allows for the occurrence of the remaining terms in Ansatz || (2). 

Replacing 8 by 8 + e where e is a small positive real constant does not change 
the spectral projection II and hence does not change (3 n . Thus 03(771) = 0. We 
postpone the evaluation of the remaining constants, co(to) and 02(771), until Section 
§ " □ 



3. Calculations on the Ball 

The eigenvalue problem on the ball is known ^] and we will only summarize 
the relevant equations for the present context. We use the following representation 
of the 7-matrices projected along ey. 



7a(m) 



\f-i ■ 7 ( m _i) 

■ 7o(m-l) 

_ y^T-l m _! 



and 



We decompose Vj = ej + 0Jj where uij = jY [ jkllk(m)ll{m) i s the connection 1 form 
of the spin connection - i.e. 

„ _ 1 (( V a \ 1 ! 

r \ \ V a / 2^ m ( m ^ a ^ m ^ 

Let P and P be the Dirac operator on the ball and the sphere respectively. In the 
notation established above, the Dirac operator on the ball is 



P 



d 771 — 1 \ If \/~^lP 

!h~ ~ 2t* / 7m(m) + r { -y/=lP 



Let (f± and Z±^ denote the eigen functions of P and P respectively, 

Pip± = ±(i<p±, PZ ( ± ] = ± + ^-j-^j Z ± ] for n = °> !» 2 > -■ 
A complete set of eigen functions is 

,„(+) c ( iJn+m^r) Z { £\SI) \ 

" i ±J n+B/ „W<»(n) J ' 

^ - a„ +rn/2 (^)^ ( o) J' where 

is the radial normalization constant. With the choice 8 = (m — l)/2 l m , the 
boundary operator A used to define spectral boundary conditions then reads 



.4 



-p q 
p 
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It is easy to determine the discrete spectral resolution of A. One can show 

z ( "\n) J V 2 / \ Z- \ty J 

for n = 0, 1, .... 



Zi n) (ft) \ f m-l \ I z m (n) 



Thus there is no zero mode spectrum. Spectral boundary conditions suppress the 
non-negative spectrum of A, which yields the implicit eigenvalue equation 

J„ +m /2-iO) = 0> n = 0,1,2, ... 

We will analyze the heat content asymptotics by considering the associated zeta 
function. Denoting by (f>k the full set of eigen functions, <pk = { l P± \ <p± ), we write 

C(*,/i,/ 2 ,AB)=X)V (h,<t>k)v (<l>H,h)&- 

fc 

The asymptotic coefficients (3 n given in equation (||) are then given by: 

(6) (3 2k (f u f 2 ,D,B) = ^(K/l/^AB), 

(7) p 2k+1 (f u h,D,B) = r(-fc-i)ResC(-A:-|,/i,/ 2 ,AB)- 

We now proceed with the explicit calculation of the heat content asymptotics on 
the ball. The ability to perform a special case calculation strongly depends on 
the choice of the initial temperature fi and of the specific heat f 2 . We establish 
Lemma ^ giving the non-locality of the heat content asymptotics by considering 
the functions: 

^^U°W md = /<2, = ( rzS o' (n) 

These spinors have the property Hf^ — whereas IT/( 2 ) = 0. Furthermore, 
because the spinor spherical harmonis involved are orthogonal, 

To evaluate the scalar product with , the relevant r-integrals are ]l4| 

f^dx x" +1 J v (fix) = ^J v+ i{p). 

We first proceed with /W. We have 

(/« ^+)) i2 =±-^ J o ^ rrm / 2 j m/2 _ l(Mr)=± I. 

Let the contour 7 enclose all positive zeroes of J m /2-i(fc)- In the contour integral 
formalism developed in [Q we have the following representation: 

C(», , , A B) = -2 / 7 f-fc- 2s - 2 |: In J m/2 _i(fc). 
Deforming the contour towards the imaginary axis we arrive at 
C( S ,/ (1) ,/ (1) ,AB) = 2/ 7 ^fc- 2 - 2 J;lnJ m/2 _ 1 (fc) 

-2^f* fc-^-^ln/^Cfc), 
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with < e G IR smaller than the first positive zero of J m /2-i(k) and with a 
semicircle around zero of radius e in the right half plane: 

7 t = {ee !t ,tG [tt/2, -tt/2]}. 

We determine the contributions of the two terms separately. Although each term 
depends on e we know the sum will be independent of e and for that reason we 
concentrate on the e-independentpart of both terms. The relevant information to 
recover the properties in (g) and (0) is encoded in the small-A: behavior of J m /i-i{k). 
The full expansion is E3 



1 OO 



from which we may obtain the expansion: 

oo 

In J v (k) = v In A; - In [2 v T{v + 1)] + ^ gi k 



1=1 



Hereby the coefficients gi are defined. In particular g\ = + 1)]. It is easy 

to see that from the small-e circle no contribution to the residues results, but that 



Ce(0, D, B) 



1. 



We consider next the contribution along the imaginary axis. This time the large-A: 
behavior is needed to determine the information in equations (|^) and (Q) . For large 
k, the relevant expansion of the Bessel function is p3 



(-1)' r(i/ + i/2 + /) 

: (2k) 1 lW(y+ 1/2-1) 



E 



and we define coefficients hj by 

1 x 
ln/„(fc) ~ k ln(27rfc) + hjk~ 



3+1 



The needed fc-integrals are trivial, dx x a = e a _ 1 , and the e-independent pieces 
at the particular values of s needed are easily obtained. We have that: 



C(o,/ (1) ,/ (1 

Res (-§,/« 

C(-i,/ (1) ,/ (1 

Res (_fe 

C(-A, A B) 



,D,B) = 0, 
,D,B) = 

,D,B) = ^ 1 (-l) fe+1 (2fc-l)/ l2fe _ 1 , fee IN, 
{-\) k 2{k-i)h 2k ^ fc-ielN. 



m — 1 
2 ' 



For the heat content coefficients we conclude from (||) and (Q) that 

/?o(/ (1) ,/ (1) ,AS) = £, /3 1 (/W,/«,AB) = -^, 
/3 2 (/W,/W,D,B) = s=±, 
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which is in agreement with 

/3 (f^J^,D,B) = J M dx(fM,f<% 

fc(fWjW,D,B) = ~f dM dyL aa {IlfW,IlfW), 

and whereby c = — 1 and c 2 = | has been determined; thus completing the proof 
of Lemma |^. As a further check that it is really the projection II that enters the 
coefficients we next consider the function f( 2 \ We have (/ ( - 2 - ) , <f±^)L 2 — ^'J 72 ^ 1 ■ 
So the starting point for the associated zeta function is 

(8) C(«, / (2) , / (2) , A B) = 2/ 7 In J m/a _x(fc). 

Note that the contour 7 encloses the zeroes of J m / 2 -i(k) only. In fact, it is possible 
to place the contour such that no zeroes of J m / 2 (k) are enclosed, because the zeroes 
of J m /2-i(k) are simple, thus ^ J' m / 2 -i( k ) = ~J m /2(k)- 
Next we use Q 

kJ m /2+i(k) = mJ m / 2 (k) - fcj m/2 _i, so 

fc2j m/2+i( fc ) = m 2 J^ /2 (k) + k 2 J^ l/2 _ 1 -2kmJ m/ 2(k)J m/ 2-i(k). 

We use the residue theorem, to see only the first term can contribute in (^). Thus 

C(s, / (2) , / (2) , A 6) = 2m 2 / 7 l^" 2 *" 4 J: In ^^(fc). 

Proceeding as before, we find 

C(0,/ (2) ,/ (2) ,fl,8) = Res C(-l/2,/( 2 )/( 2 ),D,6) = 0, 

C(-l,/ (2) ,/ (2) ,AB) = m, so 0( 2 »,/( 2 »,A6) = ^, 

/?i(/ (2) ,/ (2) ,AS)=0, and p 2 {f (2 \f (2 \D,B) = -m. 

This is consistent the form given in Ansatz [|. No boundary contributions are found 
as a result of II = 0. 
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